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ABSTRACT 

There has been a long-standing discussion in the Hterature as to whether core accretion or disk instabihty 
is the dominant mode of planet formation. Over the last decade, several lines of evidence have been pre- 
sented showing that core accretion is most likely the dominant mechanism for the close-in population of 
planets probed by radial velocity and transits. However, this does not by itself prove that core accretion is 
the dominant mode for the total planet population, since disk instability might conceivably produce and 
retain large numbers of planets in the far-out regions of the disk. If this is a relevant scenario, then the 
outer massive disks of B-stars should be among the best places for massive planets and brown dwarfs to 
form and reside. In this study, we present high-contrast imaging of 18 nearby massive stars, of which 15 
are in the B2-A0 spectral type range and provide excellent sensitivity to wide companions. By comparing 
our sensitivities to model predictions of disk instability based on physical criteria for fragmentation and 
cooling, and using Monte-Carlo simulations for orbital distributions, we find that ^85 % of such com- 
panions should have been detected in our images on average. Given this high degree of completeness, 
stringent statistical limits can be set from the null-detection result, even with the limited sample size. We 
find that <32 % of massive stars form and retain disk instability planets, brown dwarfs and very low-mass 
stars of <100 A/jup within 300 AU, at 99 % confidence. These results, combined with previous findings 
in the literature, lead to the conclusion that core accretion is likely the dominant mode of planet formation. 

Subject headings: planetary systems — brown dwarfs — stars: massive 



Introduction 



With the many hu ndreds of planet candidates from 



radial velocity (e.g. iMavor et al.l 120091: IWripht et al. 



201 lb and more t han a thousand putative candidates 



from transits (e.g. iCarbera et al.l 120091 iBorucki et al. 
201 lb . the population of close-in (~1 AU) exoplanets 



is starting to be reasonably well characterized down 
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to sub-Jovian masses. Less is known about the wide 
(>10 AU) population of planets, but the range is start- 
ing to be probed with direct imagin g. Several early 



surveys came up with null results (e.g. lLafreniere et al. 



[2007b; Kasper et al. 2007), implying that very massive 
planets and brown dwarfs are ra re in certain parameter 
ranges (" Nielsen &c"los3 12010^. However, as instru- 
ments and techniques have improved, and as a wider 
range of target classes have been examined, substellar 
companions to stars have been increas ingly discovered 
at a r a nge of orbital sep a rations (e.g. 



Th; 



20091: iBilleretaLlboiOl: lLafreniereet alJl201lh . In 



knann 



etal 



particular, some of the first probable giant plane- 
tary systems have been detected around A-stars (e.g. 



Marois et al.l 120081; lUagrange et al.ll2010l) . While the 



close-in population of planets is believed to have 
formed by core accretion (e.g. Mordasi ni et al. 2009]), 
the wide separations and high masses of some of these 
companions have led to the speculation that disk insta- 
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bilit ies may play a role in the formation of such objects 
(e.g. Bossl201 H) . The possible increase in frequency of 
massive planets around more massive stars implied by 
the direct imaging detections and shown more clearly 



in radial velocity surveys (e.g. iHekker et all |2008| 



Johnson et alj|201(} ^ may be attributed to the fact that 
more massive stars are expected to have more mas- 
sive disks, which should benefit any planet formation 
mechanism. At some stellar mass, one might however 
expect the frequency of planets formed by core accre- 
tion to turn over, because the UV radiation increases 
with stellar mass, which decreases the primordial disk 
hfetime, and which may in turn impact core accretion 
since it requires formation timescales of a similar or- 
der, a few Myrs (e.g. Alibert 2005). Disk instabilities, 
on the other hand, operate on much shorter timescales 
and should be insensitive to this effect. Hence, one 
might expect that the most massive stars should be the 
most beneficial hosts for the formation of massive gi- 
ant planets and brown dwarfs through disk instability. 

In this paper, we address the occurrence of mas- 
sive planets and brown dwarfs around massive stars 
through high-contrast imaging of 18 of the heaviest 
stars in the solar neighborhood. We examine 15 stars 
in the B2-A0 spectral type range over almost the full 
parameter range in which disk instabilities could be 
expected to form substellar companions. We investi- 
gate two white dwarfs, likely remnants of very mas- 
sive progenitors, over almost the full primordial range 
in which massive giant planets are expected to form or 
reside from core accretion. In addition, the very nearby 
A7-type star Altair is examined with a high sensitiv- 
ity to giant planet companions. In the following, we 
first present the observations and data reduction, then 
the process of data analysis, comparison with theoret- 
ical models, and Monte Carlo simulations. This is fol- 
lowed by a discussion about core accretion and disc 
instability, and the implications of this study in that 
context. Finally, we make some concluding remarks. 

2. Observations and Data Reduction 

The observations presented in this paper were 
taken with adaptive optics using primarily NIRI 



( Hodapp et al. 2003h at Gemini No rth, and addi 



tionally IRCS ("Kobay ashi et alJ l2000h at Subaru and 
NIRC2 at Keck. A summary of the observational pa- 
rameters is provided in Table [T] Our original plan for 
the NIRI data was to observe each star in both the 
CH4S- and Bra-band. The CH4S-band then provides 



optimal sensitivity at large separations where the back- 
ground noise sets the limit. At small separations. Bra 
observations would provide the best sensitivity, as the 
contrast performance in that band has been shown to 
be outstanding for bright stars w i th pre sent-day in- 



strumentation ( Janson et al. 2008, 2009). However 



during the first observations with NIRIh-ALTAIR it 
became clear that the thermal background of this in- 
strument combination is far higher than is the case 
for VLT/NACO (where the previous measurements 
had been made), hence the method could only pro- 
vide useful information for the very brightest stars in 
the sample. Instead, we used the Xc-band for small 
separations for the majority of the sample. The B2- 
AO-type stars in the sample were originally chosen 
based on the constraint that the brightness should be 
sufficient in L'-band (<4 mag) to provide an adequate 
sensitivity in the Bra-band, and based on their observ- 
ability with the given instrument in the given observing 
period. Hence, the sample is brightness-limited (and 
thus largely volume-limited for a certain stellar mass). 
Some stars were removed from the list on the basis 
of binarity that would severely affect the achievable 
contrast or dynamical stability in the systems, but no 
obvious relevant bias results from this selection. 

In all cases, the images were taken in pupil- 
stabilized mode rather than field-stabilized mode, 
to allow for high PSF stability and efficient imple- 
mentation of angular differential imaging (ADI, e.g. 
Marois et al. 2006). The PSF core was allowed to sat- 
urate for the main ADI sequences of each target. For 
the NIRI Bra- and i^c-band data and the NIRC2 H- 
band data, short sequences with the star non-saturated 
were taken adjacent to the ADI observations for photo- 
metric calibration. For the NIRI i^c-band cases where 
non-saturated images could not be acquired due to 
the brightness of the star (such as for Regulus), non- 
saturated images in the same band of other stars in 
the sample were used as photometric references. In 
the case of the IRCS data of Altair, we used non- 
saturated images of HR 8799 taken adjacent to the Al- 
tair observation as photometric standard images. For 
the NIRI CH4S-band images, a ghost feature could 
be used for p hotometric calibration as described in 
Lafreniere et al. (2007b). For each case where one or 



several companion candidates were found in the first 
epoch observation, a follow-up observation was exe- 
cuted 1-2 years afterward with the same settings (with 
two exceptions, see Sect. O. 

Reduction of the ADI data was performed with a 



2 



Table 1 : Summary of observations 



Target 


T TTFI TT~\ 

Hlr lU 


Instrument 


Filter 


Date 


Dll (s) 


Coadas 


No. rrames 


Rotation (deg) 


All- J. 

Alpheratz 


677 


NIRI 


CH4S 


2008 Oct 14 


30.0 


1 


90 


64.9 


Alpheratz 


677 


NIRI 


CH4S 


2009 Jul 30 


30.0 


1 


83 


93.0 


Alpheratz 


677 


NIRI 


Bra 


2008 Sep 22 


1.0 


25 


90 


18.3 


41 Ari 


13209 


NIRI 


Bra 


2008 Aug 29 


1.0 


25 


90 


83.5 


41 Ari 


13209 


NIRI 


CH4S 


2008 Sep 24 


30.0 


1 


90 


74.4 


41 Ari 


13209 


NIRI 


CH4S 


2010 Aug 23 


30.0 


1 


90 


64.4 


A 1 1 

Algol 


14576 


NIRI 


Bra 


2008 Sep 24 


1.0 


25 


89 


31.0 


Algol 


14576 


NIRI 


CH4S 


2008 Oct 18 


30.0 


1 


90 


34.1 


Algol 


14576 


NIRI 


CH4S 


2009 Sep 27 


30.0 


1 


45 


37.6 


A 1 1 

Algol 


14576 


NIRI 


H 


2009 Sep 27 


30.0 


1 


A C 

45 


34.2 


Bellatrix 


25336 


NIRI 


CH4S 


2008 Oct 18 


30.0 


1 


90 


46.9 


Bellatrix 


25336 


NIRI 




2008 Nov 25 


30.0 


1 


68 


33.3 


Bellatrix 


25336 


NIRI 


CH4S 


2010 Oct 07 


30.0 


1 


90 


46.9 


Elnath 


25428 


NIRI 


Kc 


2008 Nov 25 


30.0 


1 


68 


63.4 


Elnath 


25428 


NIRI 


CH4S 


2008 Nov 30 


30.0 


1 


90 


40.9 


Elnath 


25428 


NIRI 


CH4S 


2010 Oct 9 


30.0 


1 


90 


67.4 


T T Tc 1 o^r\ 

LHS 1870 


32560 


N1RC2 


H 


2009 Dec 17 


30.0 


1 


110 


58.8 


(3 CMi 


36188 


NIRI 


Kc 


2008 Nov 5 


30.0 


1 


68 


39.2 


P CMi 


36188 


NIRI 


CH4S 


2008 Nov 14 


30.0 


1 


90 


56.6 


CMi 


36188 


NIRI 


CH4S 


2010 Nov 28 


30.0 


1 


90 


43.0 


30 Mon 


41307 


NIRI 


Kc 


2008 Nov 5 


30.0 


1 


68 


27.2 


30 Mon 


41307 


NIRI 


CH4S 


2008 Nov 14 


30.0 


1 


90 


30.0 


30 Mon 


41307 


NIRI 


CH4S 


2010 Nov 28 


30.0 


1 


88 


26.2 


/ITT 

d Hya 


45336 


NIRI 


Kc 


2008 Nov 30 


30.0 


1 


68 


35.2 


6 Hya 


45336 


NIRI 


CH4S 


2008 Dec 4 


30.0 


1 


90 


40.5 


Regulus 




TVTTI>T 
IMKl 




ZUUo Dec ZU 


30.0 


1 


on 

yo 


^ 1 o 

ji.y 


Regulus 


49669 


NIRI 


Kc 


2010 Feb 23 


1.0 


30 


90 


67.4 


GD 140 


56662 


NIRC2 


H 


2009 Dec 17 


30.0 


1 


317 


91.4 


7 Crv 


59803 


NIRI 


Kc 


2009 Jan 28 


30.0 


1 


68 


18.2 


109 Vir 


72220 


NIRI 


CH4S 


2010 Jul 8 


30.0 


1 


45 


34.8 


109 Vir 


72220 


NIRI 


CH4S 


2010 Jul 8 


1.0 


30 


45 


35.3 


/3Lib 


74785 


NIRI 


CH4S 


2010 Jul 9 


30.0 


1 


45 


28.6 


^Lib 


74785 


NIRI 


CH4S 


2010 Jul 9 


1.0 


30 


45 


31.8 


eHer 


83207 


NIRI 


CH4S 


2010 Jun 29 


30.0 


1 


45 


38.9 


eHer 


83207 


NIRI 


CH4S 


2010 Jun 29 


1.0 


30 


45 


42.4 


Vega 


91262 


NIRI 


Bra 


2008 Aug 31 


1.0 


25 


189 


71.6 


A Aql 


93805 


NIRI 


CH4S 


2008 Aug 26 


30.0 


1 


68 


25.4 


A Aql 


93805 


NIRI 


Bra 


2008 Sep 06 


1.0 


25 


90 


23.0 


A Aql 


93805 


NIRI 


CH4S 


2010 Jun 16 


30.0 


1 


90 


28.1 


Altair 


97649 


IRCS 


Bro 


2009 Aug 15 


0.5 


60 


108 


29.7 
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LOCI-based ( Lafreniere et alj^ 2007ah set of routines 
in IDL, as implemented in iLafreniere et alJ (l2007bl) 
for NIRI data and with customized adaptations for the 
other cameras. Th e LOCI parameters were chosen in 
the same way as in ILafreniere et al.l (l2007bl) for all the 
data, i.e. Ns = 0.5, A^a = 300, etc. Adaptations to 
other cameras were in the form of conforming to the 
structure of the data (different keywords and conven- 
tions), but the main LOCI parameters were the same, 
although in practice there are differences due to the 
fact that the typical FWHM in units of pixels for the 
respective instruments is different (3.5 pixels for NIRI 
in CH4S and Kc, 4.0 pixels for NIRI in Bra, 4.4 pix- 
els for NIRC2 in H and 6.3 pixels for IRCS in Bra), 
which feeds into the size of the optimization regions 
and the selection of acceptable reference frames. 

All the procedures include a correction for the par- 
tial subtrac tion imposed by the LO CI algorithm, as de- 
scribed in (ILafreniere et al .l2007al) . i.e. by introducing 
several fake companions into the images and running 
the pipeline on them to evaluate the flux loss. For the 
case of 4 fim data, we intr oduced an ex t ra step anal- 
ogous to that described in ' janson et al. ( 2008h . with 



unsharp-masking with a 500 mas gaussian kernel on 
the individual frames to remove low-frequency spatial 
variations in the thermal background. 

3. Data Analysis and Results 

3.1. Identification of Companion Candidates 

Each final reduced image (see Fig. [T]for an exam- 
ple) was carefully scanned for point sources, with vi- 
sual inspection in the image itself, as well as in a cor- 
responding S/N-map where sequential annuli of the 
image are normalized by the local noise in that an- 
nulus, and where a S/N = 5 detection criterion was 
used. For each candidate, we estimated astrometry and 
photometry relative to the primary. Astrometry was 
performed using Gaussian centroiding, and photome- 
try using a circular aperture with a diameter of 77 mas, 
corresponding to a typical FWHM of the AO-corrected 
PSF in the 1 .6 /im wavelength range. The astrometric 
error was assumed to be dominated by residual dis- 
tortion correction errors as in e.g. Lafreniere et al. 
(I2OI l l). yielding errors of 6 mas in separation and 
0. 1 deg in position angle for the NIRI data. The photo- 
metric error was based on the local noise in an annulus 
at the separation of the point source. We list the prop- 
erties of each point source in Appendix A. All bright- 
ness contrast values are provided in CH4S for the NIRI 



data, except for the case of 7 Crv which was only ob- 
served in Kc- For the single point source detected in 
NIRC2 data, the contrast is in H. For this case, we 
determined astrometric errors using the pixel scale er- 
ror of 0.006 mas/pixel and the orientation uncertainty 
of 0.02° determined in previous astrometric calibra- 
tions using NIRC2 (iCihezI l2008h . No point sources 
were found in the IRCS Bra data, or indeed any of the 
Bra data from NIRI (one of the advantages of Bra is 
that background contamination is extremely low, since 
background stars are typically very faint with respect 
to the background at these wavelengths). 




Fig. 1. — Example high-contrast image of the inner 
~7"radius region of Bellatrix, with a faint point-source to 
the North in the field of view. The source is 17.6 mag fainter 
than the star. A common-proper motion test has estabhshed 
that the candidate is likely a physcially unrelated background 
star (see Fig. |2j. 

Although many of the stars had faint point sources 
in the field of view, none of these shared a common 
proper motion with the would-be host star, but were 
better consistent with the relative motion of a static 
background contaminant (with exceptions listed in the 
individual notes, including the case of Algol where 
additional observational criteria were applied for the 
companionship test). An example of a proper mo- 
tion test is shown for Bellatrix in Fig. |2] which is the 
"worst case" in the sense that it has by far the smallest 
proper motion between the two observational epochs, 
except for Algol. Two out of the three targets are fully 
consistent with the static background hypothesis. The 
remaining one is marginally deviant from the static 
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background in separation (1.9(t), but it is closer to the 
background case than to the common-proper motion 
case, from which it deviates by 2.1a. It also deviates 
from common-proper motion by 2.1 a in position an- 
gle, where it is consistent with the static background. 
We therefore consider it as a likely background star, 
possibly with a non-zero proper motion of its own. All 
astrometric points for the candidates are listed in Ap- 
pendix A. In summary, the survey revealed no substel- 
lar companions, and the data analysis was therefore 
focused on examining the implications of this absence 
of detections. 




200 400 600 800 
Time after 1st epocti (d) 



200 400 600 800 
Time after 1st epocfi (d) 




200 400 600 800 
Time after 1st epocfi (d) 



200 400 600 800 
Time after 1st epocfi (d) 



Fig. 2. — Proper motion tests for the three candidates in 
the field of Bellatrix. The upper, middle and lower rows are 
for the three different companions. The left column is for 
separation, and the right for position angle. The solid line in 
each panel is the expected motion of a background star, with 
the dashed lines corresponding to the errors. The dotted line 
marks common proper motion. The data point with errors is 
the measured position of the candidate at the second epoch 
measurement. All the candidates are probable background 
stars. 



3.2. Detection Limit Determination 

The 5a brightness contrast limits as function of sep- 
aration were evaluated by convolving each final image 
with a gaussian of the same FWHM as the observa- 
tion, and subsequently taking the standard deviation 
within a sequence of annuli as the a for each corre- 
sponding separation. In order to determine what com- 
panions cou ld have been detected, we use evolution- 
ary models jChabrier et al.lbOOOt iBaraffeet al.]|2003l) 
to translate brightnesses into masses. For this pur- 
pose, ages of the stars are required. Since B-stars 



evolve rapidly with time, isochronal dating yields their 
ages with good precision. Hence, most of t he sta rs 



were dated with theoretical isochrones ( Marigo et al. 



2008h using values of Teff and log g measured through 



spectroscopy from the literature ('Adelman & Ravlel 
2000; Gravetal. 2003; Levenhagen & Leisterl BoOet 
ISoubiran et al.H20IOl: IWu et al.ll20I ll) . In some cases 
this method was not applicable, these are noted in the 
individual notes of Sect. |3.3| for the B2-A0 stars, and 
in Sect. [33] for LHS 1870 and GD 140, and Sect. |33 
for Altair The isochrones are shown in Fig. |3] and the 
ages are listed in Table |2] In general, the targets are 
only consistent with one or two isochrones, yielding 
uncertainties of 10-20 %. Two exceptions are noted in 
the individual notes. 



Mass tracks: 2 to 9 M 




Fig. 3. — Isochrones and mass tracks from iMarigo et al.l 
in lo^TeS versus logg for the B-star spectral type 
range. Dashed lines: Mass tracks for 2-9 Mq. Solid lines: 
Isochrones for ages of 10 to 400 Myr. Also plotted are the 
error ellipses for the 14 stars in the sample that were dated 
with the isochrone method. It can be seen that from ~3 Mq 
and upwards, the isochrones are sufficiently widely spaced 
with respect to the measurement errors that isochronal dating 
can be made with good precision. 

Since the evolutionary models only provide values 
in broad-band filters by default, we used a grid of 
COND spectra dAllard et al.ll200lh for t emperatures of 
100 K to 1700 K and DUSTY spectra /Chabri er et af 



,2000 ) for 1800 K to 3000 K (surface gravities of 3.5- 
5.0 dex in both cases) to re-calculate the evolutionary 
model values from broad-band values (e.g. H) to cor- 
responding narrow-band values (e.g. CH4S) as func- 
tion of age and companion mass. 
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A common discussion regarding model inferences 
of parameters in direct imaging planet searches is 
_the discussion of hot-start (e.g. Baraffe et al. 2003t 
Burrows et al l2003h versu s cold-start (e.g. lMarlev et al. 



2007: Fortnev et all 2008) models. Cold-start mod- 



els were introduced in an effort to represent core ac- 
cretion formation of exoplanets, where the accretion 
shock may lead to signficantly cooler initial temper- 
atures than in the hot-start representation of a gravi- 
tational collapse, although there is significant uncer- 
tainty about the magnitude of this effect, due to the 
ID nature an d simplified shock tre atment of the cold- 
start models (IFortney et al.l l2008h . For most of the 
discussion of this paper, this issue is irrelevant, since 
we are considering gravitional instability as formation 
mechanism, rather than core accretion. The exception 
is in the discussion of the white dwarfs and Altair, 
where formation by core accretion may be relevant. 
However, these systems are old enough and have low 
enough minimum detectable masses that the hot- and 
cold-start models will have largely converged. Further- 
more, it is important to note that the cold-start models 
are not supported by observations, even in cases where 
core accretion seems to be the relevant mechanism. 
This is summarized in Fig. ID where temperature is 
plotted as function of age for directly imaged planets 
with constrained temperatures and small mass ratios 
relative to the par ent star, along with cold- and hot- 
start models from Fortnev et al.l (120081) . The hot-start 
models provide much better matches to the data points, 
and the cold-start models never reach the relevant tem- 
peratures, regardless of age. The most interesting case 
is that of /? Pic b, which is probably the strongest sin- 
gle case for having formed by core accretion, given 
its semi-major axis of only ~10 AU. Hence, we do 
not consider cold-start conditions as relevant for this 
study. Nonetheless, it is relevant to point out that, as 
in any direct imaging effort of very low-mass com- 
panions, the inference of physical properties from the 
observational data depends on models that have not 
been calibrated against objects with known dynami- 
cal masses, thus some caution about their accuracy is 
justified. 

The further data analysis is described separately for 
the cases of the 15 B2-A0 stars, the two white dwarfs, 
and Altair, respectively. 

3.3. The B2-A0 Sample 

The 15 stars in the B2-A0 sample are used to 
specifically address formation and in-situ retention of 



companions by disk instabilities, since the data pro- 
vide very high sensitivity in the ranges where compan- 
ions could be expected to form through such a mech- 
anism. For this purpose, we use a model that pre- 
dicts within which boundaries in mass/semi-major axis 
space fragments can form from instabilities (Klahr et 
al. 2011, in prep.). The principle behind the model 



is described in e.g. Mordasini et al. ( 201 ll) . here we 



provide a brief summary. The model requires that for 
a fragment in the planet or brown dwarf mass range to 
form, at least two conditions need to be simultaneously 
fulfilled (Rafikov 2005). One condition is that the lo- 
cal surface density in the disk needs to be high enough 
for an azimuthally symmetric instability to form at a 
certain disk radius. This can be formulated as; 



< 1 



(1) 



where Q is the Toomre parameter, Cg is the sound 
speed, K is the epicyclic frequ ency, and S is the gas 
surface density (lToomrelll98lb . In the following, this 
will be referred to as the Toomre criterion. Another 
condition is that for a fragment to form out of the in- 
stability, the cooli ng time Tcnni nee ds to be shorter than 
the shearing time (Gammie 2001), which in turn is of 
the same order as the orbital timescale Toib, i-c: 



'^cool '^orb 



(2) 



which will be referred to here as the cooling crite- 
rion. These conditions taken together mean that in a 
planetary mass versus semi-major axis diagram, there 
is a limited range where planets could form through 
gravitational instabilities. The boundaries of this area 
depend on the properties of the parent star, but as a 
general principle, the Toomre criterion states that in- 
stability occurs more easily for a higher local mass, 
so at a certain semi-major axis, only planets above a 
certain mass can form. Simultaneously, the cooling 
criterion states that fragmentation occurs more eas- 
ily for longer orbital periods, so for a given planetary 
mass, only planets beyond a certain semi-major axis 
can form. 

The model calculates the behaviour of fragmenta- 
tion based o n a ID vertical d isk structure code first 
developed in^ Bell et"aL l(ll99A but with the additional 
inclusion of realistic radia t ion inp ut from the central 
star, as in iD'Alessio et al. ( 19991) . The cooling time 
is determined by first calculating the structure for the 
isothermal equilibrium case, where the irradiation tem- 
perature Tirr, the effective temperature Toff and the 
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central temperature Tcon of the region are all equal. 
This includes a calculation of the thermal energy for 
this case, denoted £'th,o- Following this procedure, a 
second structure is calculated for a given Teen > Tcff 
with a thermal energy denoted as -Eth.i- The cooling 
time can then be calculated from a comparison of these 
cases, using the proportionality: 



ol OC 



th,l 



th,0 



crt 



(3) 



Th e most unstable wavelength (e.g. lD'Angelo et al. 
2OI0I) is used for the purpose of converting disk gas 
surface density S into the mass of the fragment. Thus, 
the model provides individual boundaries for each of 
the stars in the B2-A0 sample, taking the stellar mass, 
initial luminositjQ and metallicity as input. Metallic- 
ity is adopted from spectroscopic results in the liter- 
ature dOravet al.lliool IValdes et al.ll2004 IWu et al. 



201 Ih . In the few cases where no literature value exist 
(empty fields in Table|2]l, Solar abundance is assumed. 
A possible caveat in this procedure is if planets form 
while the star is still accreting its envelope, in which 
case the mass and luminosity would be smaller. How- 
ever, mass and luminosity have opposite effects on the 
boundaries, in the sense that a smaller luminosity low- 
ers the Toomre criterion and pushes the cooling crite- 
rion outwards, while a smaller stellar mass raises the 
Toomre criterion and pushes the cooling criterion out- 
ward. Thus, the net effect is typically modest (com- 
pa re, e.g., with the case of the Sun-like star GJ 758 
in iMordasini et al.ll201 ih . In Fig. |5]we show an ex- 
ample of boundaries for allowed planet formation by 
gravitational instability, for the case of 30 Mon. The 
mass sensitivity as a function of projected separation 
is also shown for comparison. This latter curve is a 
combination of the best detection region from each of 
the two observations performed for this target (one in 
CH4S and one in K^. To show the distribution of al- 
lowed formation regions for the full sample, we plot 
all of them in Fig. |6] Similarly, we show all the mass 
detection limits in Fig. [T] It can be readily seen that 
the detection space covers a very large fraction of the 
formation space, although since the detection space is 
presented in projected separation, the comparison does 
not provide the full story. For that purpose, we use 
Monte Carlo simulations of orbital distributions, as de- 
scribed below. 



^ In our case determined on the basis of isoclirones from lMarigo et alj 
j2008l) , at tile zero-age main sequence. 
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Fig. 4. — Temperature as a function of age for young 
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Fig. 5. — Example of model predictions and detection limit 
performance, for the case of 30 Mon. The dashed lines 
delimit the allowed range for formation, and correspond to 
the Toomre criterion and the cooling criterion, respectively. 
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cooling line are allowed. The solid line is the combined mass 
detection limit in projected separation for comparison. 
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Fig. 6. — Limits from tiie Toomre criterion and cooling 
criterion beyond tiie crossover point for all the 15 stars in 
the B2-A0 sample, to illustrate the distribution of allowed 
ranges. 
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Fig. 7. — Combined mass detection limits for all the 15 
stars in the B2-A0 sample, to illustrate the distribution of 
limits. 



In general, binaries with orbits that may cause insta- 
bility or inhibit contrast peiformance in our detectable 
image space (separations of ^0.3-10") were removed 
from our input sample. However, there are two re- 
maining targets in the sample that turn out to have 
such binary components (7 Crv and 109 Vir, see in- 
dividual notes). For these two cases, we impose the 
additional criterion that only companions with dynam- 
ically stable orbits can form in the system, where sta- 
ble orbits are determined following the procedure in 
Holman & WiegertI lll999h . Since the orbits of the two 



binaries are unknown, we assume a semi- major axis of 



1.31 ti mes the projected separation, as in lFischer et al 



(I2002h and an eccentricity of 0.36 (the m ean of the dis- 
tribution in Duquennov & Mavorlll991 ) in both cases. 
From this procedure, it follows that companions can 
only form outside of 241 AU for 7 Crv, and outside of 
102 AU for 109 Vir. 

Once the boundaries are calculated, we make two 
different lines of assumption for the distribution of 
where objects are formed and retained within these 
boundaries. One path is to assume that the distribution 
of objects is uniform over the allowed area. In terms 
of comparing to what we could detect, this is probably 
a conservative assumption, because the model refers to 
the formation of the initial fragment, which may sub- 
sequ ently accumulate fu rther gas and increase in mass 
(e.g. Kratter et al. 2010l) . which would always make it 
easier to detect in the images. We will refer to this as 
the uniform case. The other path is yet more conserva- 
tive, and assumes that the distribution is uniform along 
the minimum mass boundary set by the Toomre crite- 
rion. This coiTesponds to the extreme case in which the 
disk collapses entirely into minimum mass fragments, 
leaving no free gas to accumulate further This case 
will be referred to as the minimum case. In addition to 
the boundaries from the above description, we also set 
an upper mass limit of 100 A/jup since we are only in- 
terested in very low-mass companions (in other words, 
our study does not exclude or otherwise address bi- 
nary formation through disk instabilities, see Sect. 2). 
We also assume an outer disk radius of 300 AU. The 
outer radii of disks are typically poorly constrained, 
but millimeter interferometry observations have im- 
plied radii of a fe w to several hundr ed AU in the best 
studied cases (e.g. jPietu et alJEoOTl) . Of course, there 
is a whole range of different types of assumptions that 
could be made in addition to those we have made here. 
For instance, one might imagine that the companion 
distribution could be uniform in logarithm space rather 
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than linear space, which would give some additional 
weight to small separations which would affect the de- 
tection limit of the more distant of the targets. How- 
ever, we restrict our analysis to the uniform and mini- 
mum cases described above, as the main examples. 

The detectability of this artificial planet popula- 
tion is evaluted through use of the Monte-Carlo code 
MESS ('Bona vita etal]l201lh . The code generates 10^ 



orbits per grid point in a mass versus semi-major axis 
grid with a sampling of 5 AU in semi-major axis and 
1 Mjup in mass, for all grid points that fall into the 
allowed formation range for gravitational instability. 
The orbits are randomly oriented in space, and the 
orbital phase is also random. Since the eccentricity 
distribution is unknown, we made two separate sim- 
ulations in which the orbits were purely circular in 
one case (referred to as the circular case), and where 
the eccentricities were randomly generated with a uni- 
form distribution between 0.0 and 0.6 in the other case 
(referred to as the eccentric case). For each random 
event, the code then generates the projected separa- 
tion. The mass of the artificial planet is compared to 
the detectable mass at that projected separation from 
the mass detection limits, and the fraction of orbits at 
each grid point that turn out to be detectable then corre- 
sponds to the probability of detection at that grid point. 

As can be expected from the fact that the mass de- 
tection limits cover such a large portion of the allowed 
formation space, it is indeed the case that the proba- 
bilities of detection evaluated over the full range are 
very high. It makes virtually no difference whether 
the orbits are circular or eccentric, and the difference 
between the uniform and minimum formation cases is 
also quite small. Average detection probabilities for 
the individual stars are listed in table |3] Typical proba- 
bilities are 80-90 %. The only target with significantly 
less than 80 % consistently for all cases is Bellatrix 
(see Fig. [8] for the circular case and |9] for the eccen- 
tric case), where almost half of the artifical population 
is undetectable. The reason for this is that this star 
has the largest distance in the sample, in combination 
with the fact that it has a low metallicity, which forces 
companion formation to relatively small separations. 
For similar reasons, Bellatrix has among the highest 
differences in detection probability between the circu- 
lar and eccentric cases - this is because it is the case 
where the detection limit has the largest impact on the 
observability of the population. The reason that the 
probability is higher in the eccentric case for Bellatrix 
is that eccentric orbits with semi-major axes just be- 



low the inner working angle of the image spend most 
of their time at apastron, thus providing a net benefit in 
observability relative to the circular case. Other cases 
that have relatively big differences between the circu- 
lar and eccentric cases are Regulus and Vega, where 
the outer working angle is small enough that eccentric 
orbits at apastron near the outer edge of the field are 
preferentially lost relative to the circular case. 
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Fig. 8. — Detection probability map for Bellatrix, for the 
circular case. 
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Fig. 9. — Detection probability map for Bellatrix, for 
eccentric case. 
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Given these very high detection probabilities even 
after accounting for projection effects, some rather 
strong statistical constraints can be set on the compan- 
ion frequency from disk instability, despite the limited 



9 



Table 2: Target properties 
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size of the sample. Given that no companions were de- 
tected, it follows from binomial statistics that for a cer- 
tain intrisic frequency / of disk instability companions 
in the sample, the probability P that such a frequency 
is consistent with our observations is given by: 



(4) 



where pi is the detection probability of a given 
star in the sample. We use this procedure for all 
four combinations of cases (circular/elliptic and uni- 
form/minimum) for the 15 B2-A0 stars in the sam- 
ple. It follows that less than 29.0-31.8 % of such stars 
form planets, brown dwarfs or very low-mass stars 
of <100 Mjup within 300 AU through disk instabil- 
ity, at 99 % confidence. Here, 29.0 % corresponds to 
the circular, uniform case and 31.8 % to the eccentric, 
minimum case. Conversely, at 95 % confidence, less 
than 19.9-21.8 % form such objects and at 75 % con- 
fidence, less than 9.7-10.6 % do. 

Individual notes for some targets follow: 
41 Ari: There are no spectroscopically determined 
values for Tcff and log g for this target in the literature, 
so for the purpose of isochronal dating, we use values 
for th ese quantities determined through photometry in- 
stead dAUende-Prieto & Lamberll 1999b . 

Algol: This star is the template system for Algol- 
type interacting binaries, where mass from a secondary 
giant star overflows its Roche lobe and transfers to 
the primary. In addition, there is a tertiary (non- 
interacting) component in the system (e.g. i Kimll 19891) . 
Due to the fact that mass transfer is taking place, it 
would not be applicable to use isochronal dating of the 
primary component, since such a procedure assumes a 
constant stellar mass over time. An accurate age esti- 
mation must account for the full evolutionary history 
of the system. Hence, we use an age estimation from 
ISarna ( 1993) that does precisely this. The best-fit age 
from that procedure, which we use here, is 460 Myr 
We take the error on this quantity as the difference be- 
tween the values of the two different models calculated 
■'20 Myr. One candidate companion 



in ISarnal (119931) . 
was found in the first-epoch data of Algol. By chance, 
Algol has a very small proper motion on the sky of 
2.8 mas/yr Hence, in contrast to all the other targets, 
the motion is too small for a statistically significant 
common proper motion test over the baseline of one 
or two years. However, since the projected separation 
of the candidate is about 145 AU and the total sys- 



tem mass is ^6.2 Mq, one might expect some orbital 
motion if the candidate was real. For instance, at a 
circular face-on orbit, the motion would be 45 mas/yr 
Furthermore, if it was a real companion it would be 
~13 Afjup, and should be expected to have a signifi- 
cant non-zero value in CH4S-i?, as opposed to a typi- 
cal background star Hence, we followed up the candi- 
date in _ff-band. It shows no significant orbital motion, 
and with a contrast of 14.97±0.10 mag in H com- 
pared to 15.10±0.11 mag in CH4S, it also shows no 
dimming at all relative to the star in H. We therefore 
conclude that it is very likely a background star. 

6 Hya: While most stars can be dated isochronally 
to good precision, this is one of two targets that have 
relatively large uncertainties, since it has a relatively 
low mass and young age, where the isochrones are 
densely packed. The best-fitting isochrone is 130 Myr, 
but the fitting ages range from 60 Myr to 160 Myr 

13 Lib: Like 9 Hya, this is the other of the two 
cases where the isochronal dating has relatively large 
uncertainties, due to its relatively low mass and young 
age. The best-fitting isochrone is 80 Myr, but the fit- 
ting ages range from 40 Myr to 130 Myr 

7 Crv: The binary companion detecte d in our 



image s has been previously discovered by [Roberts 



(l2007h . but without common proper motion confir- 
mation. Comparing our data with the previous epoch, 
we can now confirm physical companionship: The 
secondary was located at 1.06±0.01"separation and 
96.8±2.0° position angle in 2002, and 1.119±0.006"and 
106.5±0. 1° in 2009. This corresponds to a motion rel- 
ative to the primary of 28 mas/yr, while the proper 
motion of the primary is 161 mas/yr, so a static back- 
ground object can be excluded with high statistical sig- 
nificance. The orbital motion is also consistent with 
the magnitude of motion one might expect for a binary 
with these properties: A system with a total mass of 
5 Mq and a 50 AU orbit has an orbital period of 158 
years. For a circular face-on orbit this would lead to a 
motion of 16° over 7 years, which is of the same order 
as the 10° motion that we measure. With contrasts 
of 7.1±0.1 mag in /-band and 5 .10±0.14 mag in K^, 
isochrones (IMarigo et al.l 120081) suggest a secondary 
mass in the range of 0.8 Mq. 

109 Vir: A likely binary companion was discov- 
ered to 109 Vir at 0.57"separation, 65.9° position an- 
gle and a CH4S brightness contrast to the primary 
of 6.04 mag. It has not been confirmed through 
proper motion tests, but due to the small separation 
and contrast, it is unlikely to be an unrelated back- 
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ground source. If i t is a real companion, isochrones 
dMarigo et al.ll2008h suggest a secondary mass of ap- 
proximately 0.65 Mq. 109 Vir is the only target in 
the B2-A0 sample with faint candidates that have not 
been followed up with proper motion tests. However, 
the two faint objects in the field would have too low 
masses to be able to form through disc instabilities (by 
a factor 6 each at their respective separations, and a 
factor 2-3 for anywhere else in the disk), so for the 
purpose of the analysis performed above, they can be 
stringently regarded as non-detections. 

Vega: Since Vega is a quite nearby system at 7.8 pc, 
the NIRI field of view only reaches out to ~90 AU. For 
this reason, we expand the field of view by including 
the contrast limits in a wider field of view provided in 



Metchev et alj (l2003h for H-band and [Macintosh et al, 



([20031) for K-band. This leads to that companions 
above approximately 7 A/jup can be excluded out to 
~230 AU. 

A Aql: This star is located only 5° from the galactic 
plane, 30° from the center, and consequently the high- 
contrast images are crowded with many background 
stars. In total 105 objects have been registered. Since 
the two epochs had somewhat different orientations 
and one had better observing conditions leading to 
a better background-limited sensitivity, not all object 
that are visible or statistically significant in one epoch 
are likewise in the other. Here we focus on astrome- 
try of those objects that are statistically significant in 
both epochs. This encompasses 55 objects, and makes 
us fully complete out to 10" and down to 17 mag con- 
trast. Such a contrast corresponds to ~4 i\/jup, which 
is more than a factor 2 lower than the lowest-mass ob- 
jects that could form through gravitational instabilities 
in the disk of A Aql. The rest of the targets are consid- 
ered unconfirmed. 

3.4. The White Dwarfs 

For the two white dwarfs LHS 1870 and GD 140, 
total ages including both estimated main-sequence 
lifetime and post-main-sequ ence cooling time were 



adopted from Burl eigh et al.l (l2002h . Primordial and 
present-day masses were also adopted from that study. 
LHS 1870 is estimated to have had a main-sequence 
mass of 3.1 A/©, a white dwarf mass of 0.66 Mq and 
an age of 700 Myr GD 140 had an original mass of 
6.3 Mq, a present-day mass of 0.9 Mq and an age 
of 250 Myr LHS 1870 has a modestly bright object 
in the field of view (see Appendix A). It has not been 



tested for proper motion, but if it were a real compan- 
ion, it would have a mass of ^15 Afjup at the system 
age, which would make it of order ^30 % as bright as 
LHS 1870 itself in M-band. This woul d make it fully 



detect able in the Spitzer excess study bv lMullallv et al 



(120071) . but their data show nothing to this effect. We 
therefore conclude that it is probably a background 
source. 

Based on our mass sensitivity curves, we perform 
Monte-Carlo simulations in the same way as for the 
B-stars for the circular case, but without imposing any 
formation limits, in order to evaluate the general detec- 
tion probabilies in mass versus semi-major axis space. 
During the mass loss phase of the white dwarfs, it can 
be expected that any planetary orbits will expand adi- 
abatically, so that the semi-major axis increases by the 
same factor that the mass decreases with. This is a 
factor 4.7 in the case of LHS 1870 and a factor 7.0 in 
the case of GD 140. In other words, a planet residing 
at 2 AU during the main sequence phase of GD 140 
would be expected to reside at 14 AU in the white 
dwarf phase. For the following analysis, we therefore 
divide the physical projected separation axis by the re- 
spective factor for the purpose of probing primordial 
semi-major axes in the system. 




4 6 
Primordial semi-major axis 



12 



Fig. 10. — Detection limits for the white dwarfs GD 140 
and LHS 1870, in the primordial separation range. 

The mass detection limits of the two white dwarf 
targets are plotted in Fig [10] It is interesting to note 
that if planetary mass scales roughly linearly with stel- 
lar mass (e.g. Alibert 2Q11J, then the mass ratio be- 
tween the star and planet may be a relevant quantity 
for what types of planets could be present in stellar 
systems of different masses, compared to our own. For 
this reason, we plot our detection probabiUty spaces in 
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primordial mass fraction versus primordial semi-major 
axis for LHS 1870 in Fig. EJand for GD 140 in Fig. 
[T2I along with the Sun-Jupiter and Sun-Saturn values. 
In this sense, a Jupiter-equivalent planet would be de- 
tectable with near-unity probability in both systems 
(98%forLHS 1870, 100%forGD 140), and a Saturn- 
equivalent planet would be detectable with 97 % prob- 
ability for GD 140. 
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Fig. 1 1 . — Detection probability map of LHS 1870, in pri- 
mordial mass fraction and semi-major axis. A mass frac- 
tion equal to the Jupiter/Sun system at the relevant separation 
would be detectable. 
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Fig. 12. — Detection probability map of GD 140, in pri- 
mordial mass fraction and semi-major axis. Mass fractions 
equal to the Jupiter/Sun and Saturn/Sun systems at the rele- 
vant separations would be detectable. 



The radial velocity study of iHekker et al.l (120081) 



has implied that massive giant planets around mas- 
sive stars (studied in their red giant phase) could be as 
much as 7 times more frequent than around Sun-like 
stars in the range of 1-2 AU (beyond this range, the 
study is incomplete). This is uncertain, since it is un- 
clear to which extent false positives contribute to the 
frequency in the former case. If taken at face value 
however, and extrapolated to 20 AU, within which 



17-20 % of Sun-like stars are estimated to host giant 
planets /Cummin g et al. 2008), it would be expected 
that massive giant planets are essentially ubiquitous 
around massive stars. Another way to approach this 
issue is to extrapolate the planet freqeuency relation in 
Johnson et al. (2010), which predicts that giant planet 
frequency is directly proportional to stellar mass, in 
which case we might expect that giant planets should 
be ubiquitous from ^5 AIq and upwards, following 
the reasoning above. In the cases of LHS 1870 and 
GD 140, we probe down to 3.5 Mjup and 2 Afjup re- 
spectively, over essentially the full primordial separa- 
tion range. Hence, down to these masses, it does not 
appear that planets are ubiquitous, although with only 
two targets studied to this high quality, it is presently 
not possible to quantify any meaningful upper limit. 

3.5. Altair 

In contrast to the B2-A0 stars, Altair is late-type 
enough that it is difficult to determine the age through 
temperature/gravity isochrones. Instead, we use the 



asteroseismological study of Altair by ISuarez et al. 



(i2005h . which yields an age range of 250-750 Myr 
We adopt 500 Myr as our baseline assumption for the 
age. 

From the mass sensitivity curve at 500 Myr, we use 
Monte Carlo simulations in the same way as for the 
white dwarfs. The resulting mass versus semi-major 
axis detection space is plotted in Fig. [T3] Also plot- 
ted are the masses and physical projected separations 
of (3 Pic b (Lagrange et al. 2010, using the mean of 



the range 8 -15 AU) and HR 8799 b-e dMarois et al 



2008, 2010*). This provides an example of the utility 
of 4 ^m data for nearby and relatively old stars - de- 
spite the fact that Altair is more than an order of mag- 
nitude older than /3 Pic and HR 8799, all the planets 
would still be detectable with this method if they had 
been present around Altair. The detection probabilities 
would be 68 % for /3 Pic b and 94 %, 96 %, 92 %, and 
76 % for HR 8799 b, c, d, and e, respectively. It is 
an interesting comparison also for the reason that the 
three stars are in the same spectral type range - Altair 
has spectral type A7V, /3 Pic A6V and HR 8799 A5V, 
according to SIMBAD. 

4. Discussion 

Since early on in the study of extrasolar planets, 
primarily two formation mechanisms hav e been dis- 
cussed in the literature: Core accretion (e.g. lPollack et al 
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19961) and disk instabilitites dBossI l2003l) . Over the 
years, evidence has been building up that core accre- 
tion is the main formation mechanism for the pop- 
ulation of close-in planets that can be studied with 
the radial velocity and transit techniques. For in- 
stance, the correlation between mejallicity and giant 
planet frequency (e.g. iFischer & Valentil l2005h . the 
fact that the giant planet radial velocity population can 
be broadly reproduced in popul ation synthesis model s 
based on core accretion (e.g. iMordasini et al. 2009), 
and the fact that the planet frequenc y increases with 
decreasing mass an d radius (e.g. Ho ward et al. 2010l: 
Borucki et al. 201 lb all point to core accretion as the 
formation scenario. This does not by itself mean that 
core accretion is the dominant mechanism for the to- 
tal planet population. Instability-formed planets are 
expected only to be able to form far outside the sep- 
aration range where radial velocity is presently sen- 
sitive. Hence, if disk instability would form and re- 
tain a larger number of planets in the far reaches of 
the system than core accretion manages closer in, 
then it would still be fair to say that disk instability 
would be the dominant mode of formation. In this 
paper, we explicitly address this question observa- 
tionally, and find that around B-stars, which should 
in principle provide among the best possible places 
for disk instability formation of massive planets and 
brown dwarfs, such objects are rare. Previous surveys 
around Sun-like stars have not quantified limits with 
regard to the particular process of instabilities, but also 
reach the general conclusion that wi de massive plan 



ets and brown dwarfs are r a re (e.g. Lafreniere et al.. 
2007bl: IChauvin et all I2OI0I iNielsen & Close! I2OIOI) . 



Meanwhile, estimations of planet frequency based on 
Kepler results in dicate a planet frequ ency of ^34 % 
out to ^0.5 AU (IBorucki et al.ll201 ll) . which is a tiny 
fraction of the separation range over which core ac- 
cretion can form planets. Given that the frequency 
of giant planets incre ases for larger separations (e.g. 
Gumming et alJ l2008h and that planet frequency in- 
creases with decreasing mass in all domains where 
it has been properly tested (e.g. Howard et al. 2010t 



Borucki et al.l 1201 Ih . it seems inevitable that core 
accretion-formed planets will be essentially ubiqui- 
tous. Put together, these pieces of evidence enable us 
at this point to conclude with some confidence that 
core accretion is indeed the dominant planet formation 
mechanism. 

This obviously does not necessarily mean that disk 
instabilities never form and retain planets. Indeed, as 



mentioned in the introduction it is the case that, al- 
though rare, massive planets do sometimes occur in 
orbits beyond where they could have been expected to 



form in-situ by core accretion (e.g. IMordasini et al. 
2009h . This does not constitute proof that they 
couldn't have formed in that way, because there 
are mechanisms by which a planet can increase its 
orbital separation, such as outward migration in a 
disk (e.g. Grida et al. 2009) or planet-planet interac- 
tions (e.g. IVeras et al.l l2009i) . but it does mean that 
disk instability should also be considered in those 
cases. Potentially interesting companions in this re- 



fe.g. Kalaset all 2005; L 


afreniei 


reetal. 2008 


) and 


above it (e.g. Billeretal.l 


2OIOI: 


Janson et al.l 


2011; 


Lafreniere et al. 2011). In manv of these cases, reg- 



ular binary formation is probably also an option - 
although high-mass ratio systems d o seem theoreti- 
cally hard to produce that way (e.g. Bate & Bonnelll 
2OO5I) . it is not yet clear whet her this is true observa - 
tionally for wide binaries (e.g. Raghavan et al. 2010l) . 
A particul arly interesting planetary case is that of 
HR 8799 dMarois et al.ll2008l l201o() . HR 8799 is in 



man y regards the most well-studied planetary system 
(e.g.ljanson et al. 2010; Hinz et al. 2010; Bowler et al. 



2010l;lMoro-Martin et al..>20ia;iHinklevi,20TTI) . but its 



formation history remains mysterious. Planets 'b' and 
possibly 'c' are beyond the range w here in-situ core 



accretion is expecte d to occur (e.g. IMordasini et al 
2009; R afikovlEoTll) . Meanwhile, planets 'd' and 'e 
are inside of the ran ge where ins tabilities could have 
occurred in-situ (e.g. lRafikovl2007.) . Hence, regardless 
of how they formed, the planets must have undergone 
migration in one direction or the other Although we 
cannot stringently exclude either scenario, we con- 
sider it of some significance that the imaged HR 8799 
system is remarkably similar to the architecture of our 
own outer Solar system, after rescaling for eq uilibrium 



temperature, as shown in lMarois et al.l (12010ft . Follow 



ing this comparison, it can be noted that Uranus and 
Neptune are also not expected to for m in-situ in the 



same core accretion-based models (IMordasini et al. 



20091) . yet few would doubt that this is the mecha- 
nism by which they formed (see the Nice model for 
a discussion on the origin of the orbits of Uranus and 
Neptune, e.g. TsiganislfeoOSh . 

A low frequency of planets formed and retained by 
disk instability does not mean that disk instabilities 
forming planetary fragments do not occur Some stud- 
ies have indicated that formed object in the gaseous 
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disk may keep growing until the y reach s tellar masses 
201(1 . 



Stamatellos & Whitworthl l2009t iKratter et al 



Also, the timescale for instabilities is short 
and may place formation in a phase where the par- 
ent star is still in a cluster environment. In such en- 
vironments, gravitational interactions between ste llar 
system s are relatively common (e.g. iMalmberg et al.l 
2007 ). Indeed, the fraction of wide binaries has been 



observationally shown to be likely affected by this 
mechanism in th e dense Orion Nebular Cluster region 
dReipurth Wide substellar companions with 

low binding energies may thus also undergo signifi- 
cant scattering Mutual interactions between multiple 
companions formed in the disk is another example of 
a mechanism that could cause scattering within the 
system or toward the field. Finally, one might imag- 
ine that planets and brown dwarfs which form far out 
in the disk migrate inwards through gaseous interac- 
tion, in the same way as expected for planets in the 
inner system. However, it would be difficult to hide 
substellar companions from detection in this manner 
Firstly, such migrati on leads to substantial g rowth of 
the migrating object. iMordasini et al.l (l2009h find that 
in the type 11 migration domain (which is what is rele- 
vant for disk instability-formed planets since they start 
out with large masses), growth occurs along lines of a 
fixed slope in log-log space, such that for every factor 
2 that the planet moves inward, it gains a factor 8.8 in 
mass. If the disk runs out of gas for the planet to grow, 
the migration stops. Hence, for instance, if a 14 A/jup 
companion forms at 50 AU around the star Alpheratz 
in our sample (the lowest mass that can form at that 
separation for this star), and migrates into 25 AU, then 
its new mass is 123 A/jup- This high mass means it 
is actually still detectable, but moreover, it is already 
well above the substellar boundary, and the limit for 
objects of interest to this study. Furthermore, even 
besides the growth issue, the rate of migration would 
have to be rather fine-tuned in order to hide the bulk of 
the gravitational instability population. Companions 
formed in the far regions of the disk have to move in- 
ward by well over a factor 5, hundreds of AU, in order 
to become potentially undetectable. Meanwhile, com- 
panions in the inner regions of the relevant space have 



^Note that if such ejections occur, the objects become iiTelevant to 
our study, since it concerns planets and brown dwarf companions, 
both of which are bound to the star by definition. For field objects in 
the corresponding mass range, such a mechanism may however be 
entirely relevant, this is not refuted or altogether addres sed by our 
study. 



to move inward by less than a factor 10, tens of AU, 
in order not to contaminate the inner planet population 
(thus diluting the brown dwarf desert and the metal- 
hcity correlation, etc.). For these reasons, we consider 
such migration unlikely. In any case, whether disk in- 
stability ever forms and retains planetary companions 
remains an open issue. 

5. Conclusions 

We have performed high-contrast imaging observa- 
tions of 18 stars, of which 15 are nearby stars in the 
B2-A0 spectral type range, 2 are white dwarfs with 
massive progenitors, and one is the A7V star Altair 
For the B2-A0 sub-sample, the resulting mass sensi- 
tivity was compared with predictions for planet and 
brown dwarf formation through disk instabilities, us- 
ing Monte Carlo simulations for orbital distributions. 
For a range of conservative assumptions of where plan- 
ets form within the boundaries allowed by disk insta- 
bility, and a range of orbital distributions, it follows 
that ^85 % of massive planets, brown dwarfs, and very 
low-mass stars below 100 Mj^p would be detectable 
within 300 AU on average. The null-detection among 
these targets therefore means that <32 % of such stars 
form and retain in-situ companions through disk insta- 
bihty at 99 % confidence, even though B-stars should 
be among the best possible places for such companions 
to form. Meanwhile, a range of results in the literature 
indicate that the frequenc y of planets formed by core 
accretion is >34 % (e.g. ' Borucki et allEoTTl) . With 
these results taken together, there is now evidence that 
core accretion is likely the dominant formation mech- 
anism for planetary companions to stars, and disk in- 
stability is not. 

The two white dwarfs GD 140 and LHS 1870 pro- 
vide a very high sensitivity to giant planets over the full 
primordial ranges where such objects could potentially 
be almost ubiquitous in the main-sequence phase. No 
companions were however detected. In the context 
of primordial star/planet mass fraction and semi-major 
axis, planets equivalent to the gas giants in our own 
system would have been detectable with high proba- 
bility in the images, had they been present. 

For the case of Altair, we demonstrated that plan- 
ets equivalent to those in the HR 8799 and /? Pic sys- 
tems would have been detectable in the system, despite 
the fact that Altair is more than an order of magnitude 
older. This is because the 4 fim technique used is well 
suited to the study of nearby and relatively old targets. 
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Eig. 13. — Detection probability map of Altair. Also plotted 
are the masses and separations of the /? Pic b and HR 8799 b- 
e planets, all of which would be detectable in the image. 
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A. Point-sources and Contrast Curves 

Here we list all point sources that were detected in the high-contrast images (Table |4|l and contrast curves for the 
various targets in different filters, at some specific separations (Table|5]l. The contrast values in Table|4]refers to CH4S, 
unless stated otherwise. 

Table 4:: Properties of all point sources. 
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z 


iz. jmu.uuD 


S9 1 1 -1-0 1 
J A. 1 1 mw. iw 


1 7 79-1-0 1 S 


901 lf\ 
ZWi w. / 


R 
D 






1 n "^8 1 -1-0 00^1 


iA(\ 'i.(\^r\ 1 

Z^O.jOzrX). IW 


1 A OA -1-0 1 9 


9008 80 
ZUUo.oU 


R 




a 


1 ^S^-l-0 OOfi 


9 /I A A'\+C\ 1 
Z'+o.^jinw. iw 


1 7S^0 1 1 


901 lf\ 
ZWi w. / 


R 
D 




1 
i 


o.jjdoihu.uuo 


^S7 08-1-0 1 
J J I .WOIlW. iw 


1 f. QO+O 1 1 


9008 Ql 


R 
D 


Z JH-ZO 


1 
i 


/;79_|_n OOfi 


^S7 0/L+O 1 
/ .WT-inw. iw 


1 7 zi7^o 1 1 


701 77 


R 


ZJH-ZO 


9 
z 


1 1 77/;_|_n OOfi 


n 78-1-0 1 
w. / omw. iw 


iO.OOIIZW. i 1 


7008 Ql 


R 


Z JH-ZO 


9 
z 


1 7 OQzL-l-O OOfi 


w.ojinw. iw 


1 7 ^Q-UO 1 1 


901 77 


R 


ZJH-ZO 




i i .uu / mu.uuu 


90S 7S+0 10 
Ayj. J jinw. iw 


1 % AQ-I-O 1 


9008 Ql 

ZWWo." i 


R 


Z JH-ZO 




1 1 808-1-0 OOfi 


9Q7 09-1-0 1 
z.y 1 .\jz.nz\j. iw 


\A 1 Q-l-0 1 


701 77 

ZWiW. / / 


R 


Z J^ZO 


/I 


1 1 70S-I-0 00^1 
11./ yj mu.uuo 


907 07-1-0 1 
Ay 1 .\J 1 mw. IW 


1J. 9Q-I-0 1 


9008 Ql 
ZWWo. V i 


R 
D 


Z J^ZO 


/I 


1 9 OOzL-l-0 OOf^ 
iz.ww^mw.wwo 


908 ^1 -1-0 10 
Ay O.J 1 mw. IW 




901 77 
ZWiW. / / 


R 
D 


Z J^ZO 


J 


89/1-1-0 OO^i 

y .oz^znw.wwu 


f,Q 70-1-0 1 
\jy . / ymw. 1 w 


1 07-1-0 1 1 


9008 Ql 
ZWWo. V i 


R 
D 


zj4Zo 


c 

J 


y.vzzitu.uuo 


f^n QA-LO 1 

/.yomu.iu 


lO.OJitU.lU 


zUlU. / / 


ID 
D 


25428 


6 


10.856±0.006 


79.54±0.10 


18.16±0.16 


2008.91 


B 


25428 


6 


10.907 ±0.006 


77.75±0.10 


18.76±0.15 


2010.77 


B 


32560 


1 


3.736±0.002 


46.40±0.02 


5.70±0.20'' 


2009.96 


U 


36188 


1 


10.472±0.006 


352.15±0.10 


14.47 ±0.10 


2008.87 


B 


36188 


1 


10.562±0.006 


352.90±0.10 


14.84±0.10 


2010.91 


B 


36188 


2 


7.541 ±0.006 


204.66±0.10 


16.55±0.11 


2008.87 


B 


36188 


2 


7.454±0.006 


204.17±0.10 


16.88±0.11 


2010.91 


B 


41307 


1 


5.820±0.006 


86.68±0.10 


15.62±0.11 


2008.87 


B 


41307 


1 


5.938±0.006 


86.17±0.10 


15.88±0.11 


2010.91 


B 


41307 


2 


9.566±0.006 


246.18±0.10 


15.60±0.11 


2008.87 


B 


41307 


2 


9.412±0.006 


246.04±0.10 


15.72±0.11 


2010.91 


B 


41307 


3 


10.080±0.006 


221.30±0.10 


16.73±0.14 


2008.87 


B 


41307 


3 


9.965±0.006 


220.72±0.10 


16.89±0.15 


2010.91 


B 


59803 


1 


1.119±0.006 


106.52±0.10 


5.10±0.14= 


2009.09 


C 


72220 


1 


5.731 ±0.006 


281.66±0.10 


16.33±0.12 


2010.52 


U 


72220 


2 


9.262±0.006 


269.45±0.10 


14.83±0.10 


2010.52 


U 



19 



Table 4:: continued. 



HEP ID 


Source no. 


Sep. (") 


P.A. (deg) 


A mag 


Epoch 


Status^ 


72220 


3 


0.572±0.006 


65.91±0.10 


6.04±0.10 


2010.52 


U 


93805 


1 


12.542±0.006 


306.40±0.10 


18.11±0.20 


2010.46 


U 


93805 


2 


12.048±0.006 


302.38±0.10 


13.27±0.10 


2008.65 


B 


93805 


2 


12.114±0.006 


302.97±0.10 


13.21±0.10 


2010.46 


B 


93805 


3 


10.784±0.006 


290.45±0.10 


17.64±0.13 


2010.46 


U 


93805 


4 


9.703±0.006 


279.84±0.10 


15.20±0.11 


2008.65 


B 


93805 


4 


9.714±0.006 


280.68±0.10 


14.92±0.10 


2010.46 


B 


93805 


5 


10.290±0.006 


270.17±0.10 


17.03±0.11 


2010.46 


U 


93805 


6 


9.484±0.006 


270.83±0.10 


18.17±0.18 


2010.46 


u 


93805 


7 


12.679±0.006 


318.19±0.10 


17.38±0.14 


2010.46 


u 


93805 


8 


10.343±0.006 


306.60±0.10 


17.88±0.14 


2010.46 


u 


93805 


9 


8.435±0.006 


283.29±0.10 


18.31±0.20 


2010.46 


u 


93805 


10 


11. 141 ±0.006 


318.11±0.10 


17.80±0.14 


2010.46 


u 


93805 


11 


9.748±0.006 


311.08±0.10 


14.20±0.10 


2008.65 


B 


93805 


11 


9.836±0.006 


3 11. 77 ±0.10 


14.02±0.10 


2010.46 


B 


93805 


12 


7.536±0.006 


285.49±0.10 


16.69±0.16 


2008.65 


B 


93805 


12 


7.575±0.006 


286.63±0.10 


16.52±0.11 


2010.46 


B 


93805 


13 


7.091±0.006 


282.88±0.10 


15.56±0.11 


2008.65 


B 


93805 


13 


7.115±0.006 


284.03±0.10 


15.57±0.10 


2010.46 


B 


93805 


14 


11.538±0.006 


324.38±0.10 


18.13±0.18 


2010.46 


U 


93805 


15 


10.414±0.006 


324.56±0.10 


15.24±0.11 


2008.65 


B 


93805 


15 


10.520±0.006 


325.08±0.10 


15.10±0.10 


2010.46 


B 


93805 


16 


9.900±0.006 


320.31 ±0.10 


17.01±0.11 


2010.46 


U 


93805 


17 


9.226±0.006 


315.22±0.10 


16.15±0.12 


2008.65 


B 


93805 


17 


9.320±0.006 


315.91±0.10 


16.12±0.11 


2010.46 


B 


93805 


18 


8.228 ±0.006 


326.77±0.10 


14.24±0.10 


2008.65 


B 


93805 


18 


8.340±0.006 


327.38±0.10 


14.32±0.10 


2010.46 


B 


93805 


19 


7.590±0.006 


320.94±0.10 


17.63±0.14 


2010.46 


U 


93805 


20 


7.845±0.006 


316.40±0.10 


18.17±0.19 


2010.46 


U 


93805 


21 


7.783±0.006 


315.21±0.10 


17.50±0.13 


2010.46 


U 


93805 


22 


6.564±0.006 


310.48 ±0.10 


17.67±0.16 


2010.46 


u 


93805 


23 


6.309±0.006 


321.23±0.10 


17.64±G.16 


2010.46 


u 


93805 


24 


5.500±0.006 


316.10±0.10 


17.18±0.14 


2010.46 


u 


93805 


25 


4.383±0.006 


274.32±0.10 


13.76±0.11 


2008.65 


B 


93805 


25 


4.384±0.006 


276.26±0.10 


13.97±0.10 


2010.46 


B 


93805 


26 


10.263 ±0.006 


260.91±0.10 


17.81±0.14 


2010.46 


U 


93805 


11 


0.95 / ztO.OUo 


OCT CO^A 1A 
ZJ j.JS±0. iO 


i /.69±0.i4 


OA 1 A A ^ 

2010.46 


T T 

u 


93805 


28 


13.412±0.006 


247.48±0.10 


16.64±0.17 


2008.65 


B 


93805 


28 


13.344±0.006 


248.02±0.10 


16.05±0.11 


2010.46 


B 


93805 


29 


13.824±0.006 


234.31±0.10 


17.86±0.20 


2010.46 


U 


93805 


30 


11.537±0.006 


236.38±0.10 


15.97±0.12 


2008.65 


B 


93805 


30 


11.449±0.006 


236.92±0.10 


15.62±0.11 


2010.46 


B 


93805 


31 


11. 103 ±0.006 


231.42±0.10 


17.25±0.12 


2010.46 


U 


93805 


32 


10.566±0.006 


227.71±0.10 


18.02±0.16 


2010.46 


U 


93805 


33 


9.238±0.006 


228.83±0.10 


18.20±0.17 


2010.46 


u 


93805 


34 


10.610±0.006 


220.22±0.10 


16.31±0.12 


2008.65 


B 


93805 


34 


10.508±0.006 


220.72±0.10 


16.13±0.11 


2010.46 


B 
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Table 4:: continued. 



HEP ID 


Source no. 


Sep. (") 


P.A. (deg) 


A mag 


Epoch 


Status^ 


93805 


35 


12.245±0.006 


213.92±0.10 


16.11±0.13 


2008.65 


B 


93805 


35 


12.127±0.006 


214.31±0.10 


16.10±0.11 


2010.46 


B 


93805 


36 


13.312±0.006 


216.58±0.10 


14.88±0.10 


2010.46 


U 


93805 


37 


8.350±0.006 


228.17±0.10 


12.51±0.10 


2008.65 


B 


93805 


37 


8.250±0.006 


228.87±0.10 


12.14±0.10 


2010.46 


B 


93805 


38 


6.115±0.006 


228.11±0.10 


17.09±0.13 


2010.46 


U 


93805 


39 


5.980±0.006 


220.52±0.10 


17.95±0.20 


2010.46 


U 


93805 


40 


5.900±0.006 


216.69±0.10 


17.40±0.15 


2010.46 


u 


93805 


41 


8.626±0.006 


206.83±0.10 


18.26±0.20 


2010.46 


u 


93805 


42 


4.422±0.006 


228.15±0.10 


16.20±0.18 


2008.65 


B 


93805 


42 


4.325 ±0.006 


229.54±0.10 


16.24±0.12 


2010.46 


B 


93805 


43 


4.277±0.006 


228.62±0.10 


16.58±0.27 


2008.65 


B 


93805 


43 


4.178±0.006 


230.02±0.10 


16.61±0.14 


2010.46 


B 


93805 


44 


8.073±0.006 


343.95±0.10 


17.45±0.13 


2010.46 


U 


93805 


45 


4.965±0.006 


337.68±0.10 


17.45±0.18 


2010.46 


U 


93805 


46 


4.730±0.006 


338.41±0.10 


17.78±0.23 


2010.46 


U 


93805 


47 


9.028±0.006 


348.91±0.10 


17.74±0.14 


2010.46 


u 


93805 


48 


9.869±0.006 


354.96±0.10 


16.74±0.14 


2008.65 


B 


93805 


48 


10.017±0.006 


355.12±0.10 


17.04±0.12 


2010.46 


B 


93805 


49 


7.644±0.006 


359.63±0.10 


18.01±0.18 


2010.46 


U 


93805 


50 


6.975 ±0.006 


358.81±0.10 


16.90±0.19 


2008.65 


B 


93805 


50 


7.122±0.006 


358.94±0.10 


17.71±0.15 


2010.46 


B 


93805 


51 


7.470±0.006 


12.42±0.10 


15.92±0.12 


2008.65 


B 


93805 


51 


7.618±0.006 


12.20±0.10 


16.19±0.11 


2010.46 


B 


93805 


52 


5.701 ±0.006 


17.43±0.10 


16.67±0.19 


2008.65 


B 


93805 


52 


5.851 ±0.006 


17.02±0.10 


17.11±0.13 


2010.46 


B 


93805 


53 


4.866±0.006 


25.29±0.10 


17.75±0.22 


2010.46 


U 


93805 


54 


9.667±0.006 


17.56±0.10 


16.27±0.12 


2008.65 


B 


93805 


54 


9.805±0.006 


17.32±0.10 


16.45±0.11 


2010.46 


B 


93805 


55 


5.248 ±0.006 


37.53±0.10 


14.82±0.11 


2008.65 


B 


93805 


55 


5.361±0.006 


36.62±0.10 


15.19±0.10 


2010.46 


B 


93805 


56 


4.841 ±0.006 


47. 47 ±0.10 


13.34±0.10 


2008.65 


B 


93805 


56 


4.937±0.006 


46.22±0.10 


13.58±0.10 


2010.46 


B 


93805 


57 


4.256±0.006 


69.57±0.10 


15.08±0.13 


2008.65 


B 


93805 


57 


4.319±0.006 


67.65±0.10 


15.43±0.11 


2010.46 


B 


93805 


58 


3.503±0.006 


92.55±0.10 


13.21±0.11 


2008.65 


B 


93805 


5o 


3. 49 J ±0.006 


AA 1 O _l_ A 1 A 

90.io±0.i0 


IT CI ^A 1 A 


OA 1 A /I ^ 

2010.40 


a 


93805 


59 


4.719±0.006 


124.92±0.10 


16.16±0.15 


2008.65 


B 


93805 


59 


4.639±0.006 


123.41±0.10 


16.31±0.12 


2010.46 


B 


93805 


60 


6.493 ±0.006 


50.51±0.10 


17.65±0.16 


2010.46 


U 


93805 


61 


6.084±0.006 


81.46±0.10 


12.14±0.10 


2008.65 


B 


93805 


61 


6.100±0.006 


80.09±0.10 


12.27±0.10 


2010.46 


B 


93805 


62 


7.030±0.006 


105.18±0.10 


15.62±0.11 


2008.65 


B 


93805 


62 


6.983±0.006 


104.07±0.10 


15.68±0.11 


2010.46 


B 


93805 


63 


5.819±0.006 


113.03±0.10 


17.46±0.15 


2010.46 


U 


93805 


64 


4.734±0.006 


204.71±0.10 


15.89±0.20 


2008.65 


B 


93805 


64 


4.612±0.006 


205.51±0.10 


16.04±0.12 


2010.46 


B 
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Table 4:: continued. 



HEP ID 


Source no. 


Sep. (") 


P.A. (deg) 


A mag 


Epoch 


Status^ 


93805 


65 


7.426±0.006 


196.74±0.10 


16.03±0.19 


2008.65 


B 


93805 


65 


7.286±0.006 


197.07±0.10 


16.14±0.11 


2010.46 


B 


93805 


66 


4.038±0.006 


170.35±0.10 


13.62±0.11 


2008.65 


B 


93805 


66 


3.899±0.006 


170.00±0.10 


13.74±0.10 


2010.46 


B 


93805 


67 


6.772±0.006 


173.61±0.10 


14.08±0.10 


2008.65 


B 


93805 


67 


6.618±0.006 


173.47±0.10 


14.03±0.10 


2010.46 


B 


93805 


68 


9.306±0.006 


172.00±0.10 


14.62±0.10 


2008.65 


B 


93805 


68 


9.171±0.006 


171.92±0.10 


14.39±0.10 


2010.46 


B 


93805 


69 


7.070±0.006 


163.80±0.10 


17.68±0.15 


2010.46 


U 


93805 


70 


6.560±0.006 


151.84±0.10 


16.31±0.14 


2008.65 


B 


93805 


70 


6.443 ±0.006 


151.34±0.10 


16.49±0.11 


2010.46 


B 


93805 


71 


10.273±0.006 


163.01±0.10 


17.51±0.13 


2010.46 


U 


93805 


72 


11.424±0.006 


165.61±0.10 


16.97±0.12 


2010.46 


U 


93805 


73 


12.287±0.006 


169.50±0.10 


16.89±0.18 


2008.65 


B 


93805 


73 


12.167±0.006 


169.45±0.10 


16.45±0.11 


2010.46 


B 


93805 


74 


12.750±0.006 


166.67±0.10 


16.23±0.13 


2008.65 


B 


93805 


74 


12.609±0.006 


166.56±0.10 


15.37±0.10 


2010.46 


B 


93805 


75 


10.116±0.006 


153.92±0.10 


17.46±0.13 


2010.46 


U 


93805 


76 


11.804±0.006 


155.30±0.10 


17.34±0.13 


2010.46 


U 


93805 


77 


12.307±0.006 


148.33±0.10 


16.73±0.11 


2010.46 


u 


93805 


78 


7.524±0.006 


115.28±0.10 


13.64±0.10 


2008.65 


B 


93805 


78 


7.456±0.006 


114.27±0.10 


13.64±0.10 


2010.46 


B 


93805 


79 


8.442±0.006 


123. 27 ±0.10 


16.34±0.13 


2008.65 


B 


93805 


79 


8.362±0.006 


122.58±0.10 


16.38±0.11 


2010.46 


B 


93805 


80 


8.481 ±0.006 


130.64±0.10 


17.50±0.13 


2010.46 


U 


93805 


81 


9.865±0.006 


135. 27 ±0.10 


16.32±0.12 


2008.65 


B 


93805 


81 


9.760±0.006 


134.75±0.10 


16.30±0.11 


2010.46 


B 


93805 


82 


7.306±0.006 


63.82±0.10 


13.77±0.10 


2008.65 


B 


93805 


82 


7.369±0.006 


62.85±0.10 


14.13±0.10 


2010.46 


B 


93805 


83 


8.060±0.006 


62.97±0.10 


15.91±0.12 


2008.65 


B 


93805 


83 


8.116±0.006 


62.13±0.10 


16.24±0.11 


2010.46 


B 


93805 


84 


7.274±0.006 


84.27±0.10 


16.47±0.14 


2008.65 


B 


93805 


84 


7.284±0.006 


83.05±0.10 


16.71±0.11 


2010.46 


B 


93805 


85 


7.829±0.006 


97.88±0.10 


16.41±0.13 


2008.65 


B 


93805 


85 


7.821 ±0.006 


96.86±0.10 


16.64±0.11 


2010.46 


B 


93805 


86 


7.157±0.006 


100.58±0.10 


17.89±0.16 


2010.46 


U 


m o AC 
93oOj 


on 


O COA_l_A AA^ 

o.j20±0.006 


TO /I O _l_ A 1 A 


1 O A /I ^A 1 A 

i/.94±U.iO 


OAAO 

2008.6 J 


a 


93805 


87 


8.562±0.006 


71.49±0.10 


13.02±0.10 


2010.46 


B 


93805 


88 


8.421 ±0.006 


78.83±0.10 


17.55±0.14 


2010.46 


U 


93805 


89 


8.288±0.006 


91.88±0.10 


18.25±0.10 


2010.46 


u 


93805 


90 


8.286±0.006 


94.21±0.10 


18.15±0.20 


2010.46 


u 


93805 


91 


9.042±0.006 


83.03±0.10 


16.01±0.12 


2008.65 


B 


93805 


91 


9.051 ±0.006 


82.16±0.10 


16.26±0.11 


2010.46 


B 


93805 


92 


10.085 ±0.006 


98.78±0.10 


14.81±0.10 


2008.65 


B 


93805 


92 


10.055±0.006 


97.98±0.10 


14.79±0.10 


2010.46 


B 


93805 


93 


11. 455 ±0.006 


84.47±0.10 


16.46±0.13 


2008.65 


B 


93805 


93 


11.455±0.006 


83.72±0.10 


16.20±0.10 


2010.46 


B 
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Table 4:: continued. 


HIP ID 


Source no. Sep. (") 


P.A. (deg) 


A mag 


Epoch 


Status^ 


93805 


94 


12.316±0.006 


67.33±0.10 


17.09±0.12 


2010.46 


U 


93805 


95 


12.014±0.006 


64.24±0.10 


14.02±0.10 


2008.65 


B 


93805 


95 


12.071±0.006 


63.63±0.10 


14.01±0.10 


2010.46 


B 


93805 


96 


13.854±0.006 


53.75±0.10 


16.63±0.16 


2008.65 


B 


93805 


96 


13. 943 ±0.006 


53.25±0.10 


16.81±0.12 


2010.46 


B 


93805 


97 


14.130±0.006 


50.72±0.10 


16.03±0.13 


2008.65 


B 


93805 


97 


14.214±0.006 


50.32±0.10 


16.34±0.11 


2010.46 


B 


93805 


98 


13.738±0.006 


47.88±0.10 


16.58±0.15 


2008.65 


B 


93805 


98 


13.836±0.006 


47.39±0.10 


16.87±0.12 


2010.46 


B 


93805 


99 


11. 381 ±0.006 


135.53±0.10 


17.26±0.18 


2008.65 


B 


93805 


99 


11.270±0.006 


135.09±0.10 


16.88±0.11 


2010.46 


B 


93805 


100 


12.784±0.006 


140.37±0.10 


17.22±0.13 


2010.46 


U 


93805 


101 


11.815±0.006 


131.91±0.10 


16.54±0.14 


2008.65 


B 


93805 


101 


11.718±0.006 


131.37±0.10 


16.44±0.11 


2010.46 


B 


93805 


102 


14.252±0.006 


132.85±0.10 


17.20±0.14 


2010.46 


U 


93805 


103 


14.679±0.006 


136.56±0.10 


17.32±0.14 


2010.46 


U 


93805 


104 


12.775 ±0.006 


348.23±0.10 


15.62±0.11 


2008.65 


B 


93805 


104 


12.919±0.006 


348.38±0.10 


14.95±0.10 


2010.46 


B 


93805 


105 


7.072±0.006 


236.74±0.10 


14.93±0.11 


2008.65 


B 


93805 


105 


6.990±0.006 


237.74±0.10 


14.84±0.10 


2010.46 


B 


"Status flag, 


where B = 


background source, C = 


confirmed companion, U = unconfirmed 







''Contrast in H-band. 
•^Contrast in ii"c-band. 
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Table 5: Contrast limits of the survey (5 a) 



T TTF» TT~\ 
HIP lU 


rllt. 


0.5 


0.15 


1 nf 
1.0 


o n/f 
2.0 


3.0 


A nff 
4.0 


5.0 


1.5 


1 r\ r\ll 

10.0 


677 


CH4S 










17.1 


18.0 


18.7 


19.6 


20.1 


677 


Bra 


7.5 


8.3 


8.7 


9.5 


9.8 


10.0 


10.1 


10.2 


10.2 


13209 


CH4S 








15.2 


16.6 


17.3 


17.9 


18.4 


18.5 


13209 


Bra 


8.4 


8.4 


8.6 


8.8 


9.0 


8.9 


9.0 


9.0 


9.0 


14576 


CH4S 










15.1 


16.2 


16.8 


18.0 


18.6 


14576 


Bra 


9.3 


9.8 


9.9 


10.2 


10.5 


10.5 


10.5 


10.5 


10.5 


25336 


CH4S 










16.6 


17.6 


18.2 


19.2 


19.6 


25336 


Kc 






13.4 


15.5 


16.6 


17.3 


18.1 


19.0 


19.4 


25428 


CH4S 










16.7 


17.5 


18.0 


19.3 


19.8 


25428 


Kc 






13.7 


15.4 


16.9 


17.3 


18.2 


19.2 


19.6 


32560 


H 


10.1 


10.9 


11.0 


11.4 


11.3 


11.3 


11.3 






36188 


CH4S 








15.0 


15.9 


16.9 


17.7 


18.7 


19.0 


36188 






11.9 


12.6 


14.8 


16.1 


16.9 


17.5 


18.4 


18.7 


41307 


CH4S 








15.0 


16.1 


16.9 


17.5 


18.0 


18.2 


41307 






11.7 


12.8 


14.6 


16.0 


16.8 


17.4 


18.0 


18.1 


45336 


CH4S 








14.0 


15.1 


16.1 


16.3 


17.2 


17.5 


45336 






11.6 


12.7 


14.7 


15.9 


16.6 


17.0 


17.5 


17.5 


49669 


CH4S 










14.7 


15.7 


16.5 


17.5 


18.0 


49669 




11.8 


12.7 


13.5 


15.5 


16.0 










56662 


H 


9.2 


9.5 


10.1 


11.0 


11.3 


11.3 


11.2 






59803 


Kc 




13.1 


13.9 


16.1 


17.3 


17.9 


18.6 


19.2 


19.4 


72220 


CH4S 


9.7 


11.3 


12.0 


14.6 


15.9 


17.0 


17.5 


18.2 


18.4 


74785 


CH4S 


10.2 


11.6 


12.4 


14.4 


15.8 


17.1 


17.7 


18.6 


19.0 


83207 


CH4S 


11.0 


12.3 


13.2 


15.2 


16.7 


17.4 


18.1 


18.8 


19.1 


91262 


Bra 


11.1 


11.6 


12.1 


12.9 


12.9 


13.0 


13.0 


13.0 


13.2 


93805 


CH4S 








15.4 


16.7 


17.4 


17.8 


18.5 


18.6 


93805 


Bra 


7.5 


7.8 


8.6 


9.4 


9.5 


9.6 


9.6 


9.8 


9.8 


97649 


Bra 


8.8 


10.5 


11.1 


13.0 


13.9 


14.4 


14.5 


14.5 


14.5 
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